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ABSTRACT 


A previously  proposed  first  order  nonlinear  differential  equation  for 
uniaxial  viscoplasticity,  which  is  nonlinear  in  stress  and  strain  but  linear  in 
stress  and  strain  rates,  is  transformed  into  an  equivalent  integral  equation. 

The  proposed  equation  employs  total  strain  only  and  is  symmetric  with  respect  to 
the  origin  and  applies  for  tension  and  compression.  The  limiting  behavior  for 
large  strains  and  large  times  for  monotonic,  creep,  and  relaxation  loading  is 
investigated  and  appropriate  limits  are  obtained.  When  the  equation  is  special- 
ized to  an  overstress  model  it  is  qualitatively  shown  to  reproduce  key  features 
of  viscoplastic  behavior.  These  include:  initial  linear  elastic  or  linear 
viscoelastic  response;  immediate  elastic  slope  for  an  instantaneous  change  in 
strain  rate;  normal  strain  rate  sensitivity  and  nonlinear  spacing  of  the 
stress-strain  curves  obtained  at  various  strain  rates;  primary  and  secondary 
creep  and  relaxation  such  that  the  creep  (relaxation)  curves  do  not  cross. 
Isochronous  creep  curves  are  also  considered.  Other  specializations  yield 
wavy  stress-strain  curves  and  inverse  strain  rate  sensitivity.  For  cyclic 
loading  the  model  must  be  modified  to  account  for  history  dependence  in  the 


sense  of  plasticity. 


Introduction 


The  subject  of  constitutive  equations  for  metals  is  presently  very  active. 
Various  approaches  are  being  pursued  which  are  intended  to  improve  on  the  body 
of  knowledge  commonly  referred  to  as  plasticity.  Although  there  are  open 
questions  regarding  the  foundation  of  plasticity  in  continuum  mechanics  theory 
the  recent  activity  was  mainly  prompted  by  discrepancies  between  the  observed 
metal  deformation  behavior  and  the  predictions  of  plasticity  theories. 

A school  of  thought  proposes  new  ideas  in  the  framework  of  classical 
plasticity  [1, 2] ; another  group  recognizes  elastic  and  inelastic  strains  which 
are  rate  dependent  [3] . Some  motivate  their  approach  by  material  science 
considerations  [4-6] . Theories  are  also  cast  in  terms  of  two  sets  of  state 
variables  evolving  according  to  separately  postulated  growth  laws.  One  set  is 
responsible  for  rheological  effects,  the  other  for  structural  effects  [7,8]. 

In  the  endochronic  theories  [9,10]  a nondecreasing  parameter  is  used  in 
an  integral  or  differential  formulation  and  it  has  been  shown  that  plasticity 
type  phenomena  can  be  reproduced  without  decomposition  of  the  strain  into 
elastic  and  plastic  parts. 

In  previous  publications  [11,12]  we  have  given  operational  definitions  for 
aging,  rate-dependence  and  history  dependence  in  the  sense  of  plasticity  and 
have  shown  that  not  all  postulated  equations  can  reproduce  what  are  commonly 
considered  plasticity  effects.  We  have  further  argued  that  plasticity  effects 
are  observed  due  to  the  internal  changes  in  the  material  affected  by  the  loading. 
Accordingly  we  postulated  that  a measure  of  internal  change  'separate  from 
repositories  for  rate  dependence  and  aging)  should  be  introduced  in  a consti- 
tutive equation  for  metals.  In  [13]  the  properties  of  a continuously  varying 
structure  parameter  together  with  an  integral  constitutive  equation  were  in- 
vestigated and  it  was  shown  that  typical  plasticity  responses  (cross-hardening, 


cyclic  hardening,  Bauschinger  effect)  could  be  modelled.  However,  it  was 
necessary  to  solve  an  integral  equation  to  obtain  the  growth  for  the  struc- 


2. 


ture  parameter. 

A discontinuous  growth  law  for  the  structure  parameters  was  proposed 
in  [11].  It  was  argued  that  in  a continuum  approach  plasticity  effects  can 
only  be  recognized  after  unloading  and  reloading.  According  to  this  view 
plasticity  and  viscoplasticity  are  indistinguishable  from  nonlinear  elasticity 
and  nonlinear  viscoelasticity,  respectively  as  long  as  we  consider  only  loading 
(see  [11]  footnote  on  p.75). 

Here  a previously  proposed  differential  constitutive  equation  is  repre- 
sented as  an  equivalent  integral  equation.  The  integral  formulation  permits 
the  study  of  the  limiting  behavior  under  large  strains  and  large  times.  We 
demonstrate  by  qualitative  arguments  and  by  numerical  experiments  that  the 
proposed  "equation  of  state"  can  represent  many  key  features  of  viscoplastic 
behavior  as  long  as  there  is  no  substantial  reversed  loading  involved.  For 
reversed  and  cyclic  loading  the  proposed  nonlinear  viscoelasticity  law  must 
be  extended  by  the  introduction  of  new  origins  and  by  updating  the  material 
parameters  [11, 14] . A three-dimensional  formulation  of  the  present  theory 
is  given  separately. 

The  Proposed  Constitutive  Equation 

For  the  uniaxial  case  we  previously  proposed  [14, 15) 

m[a,e]e  + g[e]  = a + k[a,e]a  (1) 

as  the  simplest  constitutive  equation  capable  of  reproducing  rate  sensitivity, 
creep  and  relaxation.  In  the  above,  ct  denotes  the  axial  component  of  the 
Cauchy  stress  (true  stress),  e is  the  engineering  strain, and  e + 1 is  the  axial 
component  of  the  deformation  gradient.  The  square  brackets  following  a 
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symbol  denote  "function  of"  and  the  dot  over  a symbol  denotes  the  material  time 
derivative.  Equation  (1)  is  valid  for  infinitesimal  and  finite  strain.  The 
functions  m[  ] and  k[  ] are  positive , bounded  and  are  even  functions  for 
both  variables  taken  together,  i.e. 


and 


m[a,  e]  = m[-a,  -e] 
k [o',  e]  = k[-a,  -e]  . 


} 


(2) 


The  function  g [ ] is  odd.  The  functions  m [ ],  k[  ] and  g[  ] are  assumed  to 
have  continuous  derivatives. 

The  above  stated  requirements  are  derived  from  considering  the  qualita- 
tive behavior  of  materials  in  tension  and  compression  and  the  responses  pre- 
dicted by  (1).  For  the  actual  description  of  a particular  material,  the  three 
material  functions  m[  ],  k[  ],  g[  ] must  be  determined  from  experiments  [16]. 

Integral  Formulation 

The  constitutive  equation  (1)  can  be  expressed  as  the  following  equiva- 
lent integral  equations 

°[T1  ■ 9Ie[T”  + J - s'UMi} 

O 

dx 


exp 


(-J 


k [o[x] , e [x]  ] 


) 


s[t]  dT 


(3) 


« [tl  - J + SI 11}  exP(-  f — „ ds)  dr  (4) 

L m[  ] m [ ] J J e [s]m[o[s] , e [s]  ] / 

o T 

where  we  have  assumed  o(t  = 0)  =0  and  e(t=0)  =0.  In  the  above,  ' denotes  the 

. * 
derivative  of  a function  with  regard  to  its  argument  . 


In  the  case  of  discontinuous  strain  (stress)  rates  we  propose  (3)  and  (4)  as 
fundamental  forms. 


4. 


It  does  not  appear  to  be  possible  to  represent  (1)  by  an  equation  of  the 
type  t 


ct[t]  = J G [e  [t] , t - t]  e [t]  dT 
o 

t 

e [t]  = J J[CT[T],  t-T]&[T]  dT  . 


(5) 


(6) 


Equations  (3)  and  (4)  simplify  to  the  regular  convolution  integrals  of  linear 
viscoelasticity  if  g[e]  is  linear  in  e and  m[  ] and  k[  ] are  constants. 

Behavior  in  Tensile  Tests 


We  now  use  (3)  to  simulate  a tensile  test  under  constant  strain  rate 
e = a such  that 

e = at  . (7) 

We  designate  the  response,  i.e.,  the  solution  of  (3)  subject  to  (7)  by 


q[e]  = ate  ;<*]  = cr[t] 


e=at 


(8) 


and  a is  considered  a parameter. 


Substitution  of  (7)  into  (3)  using  (8)  yields 

e e 

a[e  ;ar]  = a[e]  =g[e]  + J - g'tx]  lexp(-  ± [ 

o k[a[x],x]  01  x 

For  extremely  slow  loading  a -*  0;  we  obtain  from  (9) 

5[e;o]  = lim  ate;®]  = g[e] 

Qf-0 

and  for  extremely  fast  loading  a — °°,  (9)  furnishes 


ds 


k [a[s] , s] 


dx  . 


(9) 


(10) 


tcr[x],x] 


(ID 


Appendix  I shows  the  parallel  development  for  (4) 


5. 


From  (10)  we  see  that  g[e]  furnishes  the  stress-strain  curve  which  is 

obtained  for  extremely  slow  loading.  We  designate  the  stress-strain  relation- 

* 

ship  for  infinitely  fast  loading  by  g[e],  From  (11)  we  have  for  infinitely 
fast  loading 


[£]]  = 2'[e] 
k[o[e],[el] 


(12) 


The  relation  (12)  between  m[  ] and  k[  ] is  only  required  for  fast  loading. 

However,  as  a useful  and  physically  motivated  simplification  of  the  theory, 

we  propose  (12)  to  be  true  for  all  loadings.  Consequently  (9)  is  rewritten  as 

e e 

a [e ; o’]  = ct [e ] = g [e ] + f jg  ' [x]  - g'  [x] }exp(-  — f — ) dx  . (13 ) 

& L J v a " k [o[s] , s]' 


Using  (12)  we  can  also  rewrite  (1),  (3)  and  (4)  accordingly. 

From  (13)  we  see  immediately  that  for  every  e 2;  0 and  a ^ 0, 


g [ e ] s ofe.-Qf]  ^ g[e] 


(14) 


provided  g [e]  - g [e]  2 0.  Therefore,  the  response  for  0 s a s <®  is  between 
the  extremely  fast  and  the  extremely  slow  response.  Moreover  for  a?  > > 0 

we  have  for  every  positive  e 


ct  [e  ;o-2 ] s CTfe,-^] 


(15) 


so  that  the  stress-strain  curve  for  the  fast  rate  is  above  the  curve  for  the 
slow  rate,  as  shown  in  Fig.l.  The  above  behavior  reflected  in  the  relations 
(14)  and  (15)  would  be  reversed  by  setting  g'fe]  - g'[e]  so  for  e s 0. 

For  or  s 0 the  responses  are  symmetric  with  respect  to  the  origin  by 
virtue  of  (2)  and  the  fact  that  g is  odd. 

Differentiation  of  (13)  with  respect  to  e yields 

e e 


do 

de 


= g [e] 


1 

ork[a[e],e] 


Jfa ' tx]  - g ' [x] }exp(-  ^ J 


ds 

k[o[s],s] 


(16) 


6. 


Performing  the  limits  for  O'  -*  0 and  O'  -*  <*>  in  (16 ) we  obtain  results  consistent 
with  (10)  and  (11),  (12),  respectively.  Moreover  we  see  that 

t*  = 01  = g'tO]  (17) 

for  all  a.  All  the  curves  in  Fig.l  have  therefore  one  initial  common  slope 

* 

equal  to  the  initial  slope  of  g irrespective  of  the  loading  rate. 

Behavior  at  Large  Strains 

We  are  interested  in  the  behavior  of  our  model  at  large  strains  in  a 
tensile  test.  To  this  end  we  investigate  (16)  for  the  mathematical  idealiza- 
tion e oo  for  a = const.  . The  resulting  indeterminate  expression  is 
resolved  (see  Appendix  II)  to  obtain 

A 

da  , 

lim  = 9 M (18) 

e — »co  at 

for  all  finite  a and  provided  g'M,  g'[°°]  and  exist  and  are  finite  . 

This  is  a very  strong  result  requiring  that  all  response  curves  become  parallel 
to  each  other  for  large  strains.  Figure  2,  which  represents  the  results  of 
the  numerical  integration  of  (13),  demonstrates  this  behavior.  For  convenience 

k was  selected  to  be  constant.  Initially  linear  elastic  behavior  is  obtained 

/ */  * 
as  g [0]  = g [0] . The  representation  of  the  g and  g functions  follows  the 

development  of  the  y(x)  functions  using  the  second  kernel  of  [22], 


While  these  mathematical  limits  may  seem  unrealistic  they  are  indeed  useful 
because  the  solutions  of  the  equations  are  rapidly  asymptotic  to  these  limits; 
see  Fig. 2. 

Throughout  this  paper  we  consider  only  k- functions  which  are  bounded.  In 
the  case  where  k[a,  e]  =k[cr-f(e)]  for  some  function  f,  we  similarly  require  that 
k[<=]  be  finite. 
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Performing  the  limit  for  (13)  in  a similar  manner  we  obtain 

lim  (o-g[e])  = (g'[®l  - g ' [°°]  )ak  [a[°°]  ,<*>]  . (19) 

e-*00 

The  limiting  value  of  (a-g[e})  is  finite  and  is  determined  by  the  strain  rate 

/ * / 

and  the  values  of  the  functions  k [ ] , g [ ] , g [ ] at  infinity. 

The  relations  (18)  and  (19)  were  derived  for  constant  strain  rate  loading.  j 

By  starting  with  (3)  and  investigating  the  limit  for  t -*  ® the  following 

1 

expressions  result: 

lim  o[t]  = g ' [e  [°°]  ] e [®]  (20) 

t-“  ' 

and 

lim  (C7  — gfeft]])  = (g ' [e  [“]  ] - g ' [e  [«]  ] ) e [®]  k [a  [«]  e [»]  ] (21) 

t-® 

provided  s[»]  together  with  the  limits  of  the  other  functions  exist  and  are 
finite. 

If  e [oo]  = “ the  expressions  (20)  and  (21)  are  identical  to  (18)  and  (19), 

respectively,  the  only  exception  being  that  a[e],  the  response  to  constant 

strain  rate  forcing,  appears  in  (18)  and  (19)  whereas  the  response  in  (20) 

and  (21)  is  a[t] . Therefore  the  responses  for  large  strain  and  time  under 

variable  e become  identical  to  the  ones  under  constant  e = o'  provided  e [™]  = o' 

and  it  is  not  necessary  to  have  always  constant  strain  rate  on  0 $ t ^ ® to 
da  / 

arrive  at  lim  — = g [e  [«>]]. 

Instantaneous  Change  of  Strain  Rate 

We  consider  two  tensile  tests  with  respective  strain  rates  a 2 and  with 

A A 

ck^>  0i^>  0.  Then  from  previous  arguments  we  know  that  ^ a[e,a^]. 

Suppose  instead  that  we  initially  load  with  and  at  some  time  tQ  and  associ- 
ated strain  we  instantaneously  change  the  strain  rate  from  to  as 
shown  in  Fig. 3.  We  call  ^[6]  the  response  to  a 2 and  &B[e]  the  response  to 
or^  followed  by  0'2  as  shown  in  Fig.  3. 

I 


8. 


First  we  investigate  the  cj  vs.  e relationship.  For  e £ e , crA  [ e]  s crn[e]  by 


previous  arguments.  For  e s the  reo^ective  responses  are 


and 


t 6 

5ft[e]=gte]  +J  (g ' [x]  - g ' fx]  )exp(- 


ds 


) 


2 x k[aA[s],s] 

e e 

a fe]=g[e]  + [*  (g  ' [x]  - g'  [x]  )exp(-  — [ — J 

B J v J k[d  [S],s]/ 


dx 


dx 


+ (ct  [g  ] - g [e  ] )exp( — — [*  — — — ) . 

Bo  o ' p k [cr^  [s] , s] ' 


e ''  b 
o 


(22) 


(23) 


where  for  e £ e a [e]  = a [at  - (a?  -ar_  )t  ] = cj[t]. 

O D D Z 4 I O 

For  large  strains  (e  -*  ®)we  obtain  from  (22) 

lim  (ct  - g [e] ) = {g'[°°]  - g ' [“]  }ar  k [5  [°°] ,”] 

g— *00  A Z Pi 


and  from  (23) 


lim  (ctr  - g [e] ) = {g'M  - g ' [°°]  }a2k  [&b  [«]  j00]  . 


(24) 


(25) 


Unless  k[cr  C00 ] , 00 ] = k[o _ [°°] , °°]  the  two  paths  on  the  ct  - e plane  will  not  coin- 
cide  for  large  values  of  strain.  In  both  cases  the  slopes  approach  g ' [<»]  as 
e -*  <*>.  We  see  therefore  that  the  two  curves,  i.e.,  ct,  and  ct  , become  parallel 
at  large  strains. 

The  question  arises  whether  the  path  ct  can  be  crossed  by  path  ct  . From 

A B 


(1),  (12)  and  the  chain  rule  we  obtain 


dCT  */  r i 

dl  = g [S]  " 


a-  g[e] 
k [CT  [G  ] , £ ] O' 


(26) 


From  Eq.  (26)  we  see  that  for  a given  stress,  strain,  and  strain  rate, 

dCT 

there  exists  a unique  slope  — irrespective  of  the  prior  history.  Since  cross- 


ing of  the  ct  and  ct  curves  would  require  a nonunique  slope  at  the  crossing 
A B 


9. 


point,  the  condition  necessary  for  crossing  cannot  be  met  and  ct  s ct  always; 

A 6 

the  two  curves  can  never  cross*. 

To  investigate  the  behavior  at  large  times' we  start  from  (3)  and  obtain 
(note  ct  and  a are  the  responses  in  cr-time  coordinates) 

A B 


lim  (ct  - g[a  t] ) = (g'M  - g'  [=°]  ) a_k [ct  [®] 
t— » oo  a z z Pi 


(27) 


and 


lim  (ct  -g[«t  + Of  (t  - t ) ] ) = (g / [“>]  - g 7 [“])  a_k  [CT  [®  ] , °°]  . (28) 

■£- — » oo  d ± O Z O ZB 

If  g 7 [ro]  = E then  we  have 

b 

lim  (ct  - ct  +E  t (or  - or  ) ) = (g7  [°°]  - E )a  (k[cj  [°°] ,“]  - k[cr_  [®] ,°°] ) . (29) 

^_*co  A B o O Z ± b Z Pi  B 

Even  if  the  right-hand  side  of  (29)  would  be  zero,  ct  [»]  - CT  [<*>]  # 0 i.e.,  the 

A B 

graphs  of  ct  vs.  time  and  cr  vs.  time  do  not  coincide.  The  corresponding 

CT-e  graphs,  however,  may  eventually  coincide  as  can  be  seen  from  (24)  and  (25). 

So  far  we  have  only  considered  > 0.  The  previous  arguments  can 

be  extended  to  the  case  of  a,  > a„  > 0 with  the  result  that  the  curve  for  a 

12  B 

vs.  e will  approach  ct  vs.  e from  above  but  will  never  cross  the  ct  curve. 

A A 

Change  of  Slope 

The  slopes  at  point  A in  Fig. 3 approached  from  the  left  and  from  the 
right  are  determined  from  (16)  or  (26)  to  be 


..  ^_*,r  . %[eo]-g[eo] 

e-iV  de  9 e°  a k[CT  [e  ],e  ] 
o 1 B o ’ o 


and 


d°B  *, 

lim+ 


e—  e 


Veo]  - g[eo] 
a2kl5B[eo]’eo] 


(30) 


(31) 


This  argument  is  due  to  Dr.  H.  Moon. 


10. 


so  that  the  slopes  differ  by 


do  do 

*snr  <e  = 0 --5T  (e  = V 


a2-ai  aB[eo]-greo] 


2 1 ' ktaB[eo],  eQ] 


(32) 


We  obtain  therefore  an  instantaneous  change  in  the  slope  of  the  a - e curve  as 

the  strain  rate  is  changed  instantaneously  from  a to  or  . For  a > or  the 

dcr  1 do  2 

T>  Jf.  Q «» 

right-hand  side  of  (32)  is  positive  so  that  — — (e  = e ) > — — (e  = e ) . For 

as  o de  o 

a;^  > > 0 the  opposite  is  true. 


Behavior  at  Small  Stresses  (Strains) 

Here  we  are  interested  in  the  behavior  around  the  origin.  Suppose  that 

i * 

for  e £ e (e  >0)  the  functions  g and  g can  be  approximated  by  g[e]  = Ee 

3 3 

* * 

(E  is  Young's  Modulus)  and  g[e]  = Ee  and  that  m[  ] and  k [ ] are  constant  for 
all  | e | <■  e and  associated  |a[  £ a (a  >0)  and  equal  to  m and  k , respect- 

3 cl  3.  3 a 

ively.  Then  from  (3)  for  |e|  ^ e 

3 

t 

a[t]  = Ee  + J (E  - E)exp(-  ~ (t-T)J  e dT  . (33) 

o 3 

* * 

For  E > E (33)  is  the  standard  linear  solid  of  viscoelasticity.  For  E = E 

initial  linear  elastic  behavior  is  reproduced  on  |e|  £ e^,  even  for  nonconstant 
m and  k.  These  two  models  are  important  representations  of  solids  and  can  be 
easily  obtained  from  (1)  just  by  constructing  the  coefficient  functions 
appropriately  [22] . 


Creep  and  Relaxation 

We  use  Eq.  (3)  to  compute  the  stress  response  in  a relaxation  test  where 
the  strain  is  increased  to  some  value  e and  then  kept  constant  such  that 
e = 0 . We  assume  the  strain  history  to  be  given  by 


11. 


e = e [t] 
e = 0 


0 £ t < t 


t ^ t 


with  e [t  ] = e 
o o 


where  e [t]  is  some  smooth  function  of  time.  From  (3),  using  (12),  we  obtain 

for  the  relaxation  stress  ct  for  t ^ t 

R o 

t 


Vt]  =g[eo]  +(exP  - J k [o  ‘[x] , e ] ) <Veo>  - g[eo]) 


where 


t ' R 
o 


l 

aR[eo]  -g[eo]  = J Cg/[efT)l-g/r®(T)]}exp{  J ■ ^ [x]]}e 

*-  R ’ 


►e  (T)dT  . 
(34) 


Alternatively  if  e [t]  = e h[t-t  ] and  e[t]  = e 6[t-t  ] where  h[t]  and 

o o o o 

6[t]  are  the  step  function  and  the  Delta  function,  respectively,  i.e.,  if  we 


consider  an  instantaneous  loading  from  zero  to  e at  time  t we  obtain  from  (3), 

o o 


for  t S t , respectively, 


t 

r* 


CTfi[t]  =g[eQ]  + {g[eol  -g[eQ]}  exp  - J 


dx 


J k [a. [x]  e ] 
t 0 7 o 


(35) 


where  [t]  is  the  relaxation  stress  due  to  the  impulsive  strain  input. 

Differentiation  of  (34)  and  (35)  shows  that  the  relaxation  rate  is  always 

*/  / 

negative  or  zero  provided  g -g  ^0.  For  large  times  t ~*  <*>,  a[t-*°°]  = g Te Q] 

for  both  (34)  and  (35).  The  curve  g[e]  is  therefore  not  only  the  slow  loading 

curve  but  also  the  terminal  curve  of  all  relaxation  tests. 

From  (35)  we  see  that  the  impulsive  stress  relaxation  starts  at  t=t 
* 

from  g[eQ],  i.e.,  from  the  curve  traced  out  in  the  infinitely  fast  constant 
strain  rate  loading.  If  the  loading  up  to  eQ  is  smooth,  such  as  by  constant 
strain  rate,  then  we  observe  from  (34)  that  the  relaxation  starts  from  a point 
on  the  curve  traced  out  in  the  constant  strain  rate  test. 
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The  decrease  in  stress  during  the  relaxation  test  is  controlled  by  the 
function  Jc  [ - ] . We  may  therefore  call  it  the  relaxation  function.  To  demon- 
strate the  effect  of  the  function  k on  relaxation  behavior  we  have  computed 

* 

the  relaxation  from  a point  on  g in  Fig. 4.  Curve  1 corresponds  to  con- 
stant k and  (34),  (35)  are  linear  in  stress.  For  curve  2 k = k[CT-g[e]]  as  in 

Fig. 10  so  that  (34)  and  (35)  are  nonlinear  in  stress.  In  Fig. 4 the  time 
scale  depends  uf)on  the  time  scale  used  for  k;  if  k is  given  in  units  of  hours 

then  the  abscissa  is  measured  in  terms  of  hours. 

Appendix  III  gives  a discussion  of  the  creep  equation  obtained  from  (4). 

A comparison  of  (34)  with  the  corresponding  creep  equation  (All)  shows  that 

the  exponential  decay  in  a creep  test  is  controlled  by  m[  ],  g[  ] and  e[  ] 
and  involves  a complicated  interaction  of  these  quantities.  Although  creep 
and  relaxation  show  an  exponential  type  of  decay,  their  detailed  evolution 
is  quite  different. 


Spacing  of  Creep  and  Relaxation  Curves 

Consider  a series  of  creep  (relaxation)  tests  performed  at  an  increasing 

sequence  of  constant  stress  (strain)  values;  we  require  that  the  resulting 

curves  of  strain  vs.  time  (stress  vs.  time)  and  of  strain  rate  vs.  time 

(stress  rate  vs.  time)  are  nonintersecting  as  is  usually  the  case  with  actual 

test  results.  Further  we  require  that  any  two  curves  e 2 [t]  and  e^tt],  (ct0 [t] 

2 1 

and  [t]  ) associated  with  two  constant  stresses  cj  > a (constant  strains 
1 o o 

2 1 

e > e ) satisfy  the  condition  e „[t]  s e,  [t]  for  any  time  t (a_rt]  s cr,  [t] ) 
o o 2 1 2 1 

and  that  e2[t]  s fe^tt]  (&2 [t]  2 a [t]).  These  conditions  require  that 


Se[t]  ^ „ , Sett] 

— r ^ o and  — 

Scr  S O 


> 0 


(36) 
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a 0 and  a 0 (37) 

S€o 

for  the  creep  and  relaxation  test,  respectively;  eft],  eft]  in  (36)  depend 
parametrically  upon  ct^,  and  aft],  ft]  in  (37)  depend  parametrically  upon  €q. 

If  we  specialize  (1)  for  the  creep  test  (a  = 0,  a = a ) and  perform  the 
differentiation,  with  respect  to  c , we  obtain 

Fit]  - Hit] 

: — rb — 

o m[CT  , e ft]  ] 

o7 


where 


F ft]  = m [ctq,  eft]]  - (ctq  - g [e  ft]  ] ) 


dm  [a  e ft]  ] 

O' 


dm  [a  , eft]] 

H ft]  = m [a  e ft]  ]g'  fe  ft]  ] + (a  - g [e  ft]  ] ) 


Equation  (38)  is  integrated  to  obtain 


de  ft]  = f F [t  ] 

^CT  J far  r -.I2 

° O m[CT  e[T]] 


(exp  - % 


dx  dT 


t m[CT  , e [x]  ] 
o7 


+ aap  . J 

o O info  e [x]  ] 


where  in  (41)  t = 0 corresponds  to  the  initiation  of  creep,  e [0]  # 0. 
de  ft  = 0] 

When  — r 2 0 from  (38)  and  (39)  we  see  that  sufficient  conditions  for 

0(7  — 

O 

satisfying  (36)  are 

F ft]  2 0 and  Fft]  - H ft]  2 0.  (42) 

o 

Following  similar  procedures  for  relaxation,  sufficient  conditions  for 


satisfying  (37)  are,  respectively 
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dk[<7[tj,e  ] 

P[t]  =k[CT[t],eo]g/[eol  + (a[tl  -g[e  1)  - 


9e 


s 0 


and 


with 


Ptt]  " Q[t]  a ° 


Q[t]  = k[a[t] , e ] - (a[t)  -g[e  ]) 
o o 


9k  fa  [t]  e 1 
o 

9a 


(43) 


(44) 


(45) 


The  examination  of  (39)  - (42)  shows  that  the  no  crossing  conditions  (36)  and  (37) 
are  in  general  not  satisfied. 

From  Eqs.  (38)  - (41)  and  their  counterparts  for  relaxation,  we  obtain  the 
limits  for  large  times 


. . 9e  [t]  _ . ( F[t] 

ii®  ti®  TlTi^i] 


) _ i 

/ g ' [e  [a 


lim 


9a  ft] 


t“*oo  9c 


lim  = 0 

t-®  Bcto 


■ (frlf) 


lim  = 0 . 

t-®  de0 


(46) 


(47) 


(48) 


(49) 


where  we  have  assumed  that  in  creep  or  relaxation  (a-g)  -•  0 as  t -•  ®.  There- 
fore (46)  and  (47)  do  not  hold  for  those  creep  conditions  for  which  there  is 
no  solution  to  the  equation 


a = g [c] 
o 


(50) 


i.e.,  in  those  cases  for  which  (1)  reproduces  secondary  or  tertiary  creep; 


see  also  Appendix  III. 
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Useful  Specializations  of  (1) 

A)  The  Functions  m [ ] and  k[  ] are  Selected  to  Depend  Upon  (a-g[e]) 

*,  $ 

and  g [e]  =E  (Modulus  of  Elasticity) 

Because  of  this  special  dependence  a and  e have  no  independent  influence 

on  the  behavior  of  (1).  For  a given  loading  rate  the  properties  of  (1)  are 

determined  by  the  value  of  CT-g[e]  alone,  and  curves  equidistant  to  a=g[e] 

are  curves  alorig  which  the  properties  of  solutions  of  (1)  are  identical. 

Specifically  with  reference  to  (26)  we  see  that  the  slope  ^ for  a given 

constant  e is  constant  on  a curve  a-g[e]  = const.  Moreover  for  a-g[e]  =0 
da 

— = E irrespective  of  the  strain  rate.  The  slopes  of  all  of  the  curves  de- 
parting from  the  curve  a = g[e]  are  all  equal  to  the  elastic  modulus.  (It  is 
however,  not  possible  to  depart  from  a-g[e]  =0  by  creep  or  relaxation,  as 
can  be  easily  verified  from  (1).)  With  the  a-g[e]  dependence  in  the  functions 
the  curve  a = g[e]  assumes  special  significance  since  along  it  the  behavior 
predicted  by  (1)  is  the  same  as  at  the  origin. 

While  every  solution  curve  of  (1)  for  constant  or  variable  strain  rate 
must  become  parallel  to  g[e]  [see  (18)  and  (20)]  only  the  a-g[e]  dependence 
in  k models  a nonlinear  dependence  between  a and  a-g[e].  Consider  now  two 
tests  with  and  O' 2 and  associated  stress  responses  a^  and  a 2,  respectively. 
From  (19)  we  obtain  the  ratio 

(a^”]  - gM  )k[a2  M - gM  ] 


(&2  M - g M ) k [a^  M - g M ] °2 


and  therefore  (CT^-gM)  and  (c^-gt00])  are  related  in  a nonlinear  fashion. 
For  any  other  choice  of  the  a,  e dependence  in  the  functions  m [ ] and  k [ ] 


We  recall  the  conditions  k(°°)  = finite  and  k positive. 
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we  have  for  [m]  = [»]  = ®,  k [6^  [<*>] , ®]  = k[&2  [“] , °°]  *so  that  the  ratio  (51) 

reduces  to 

a [“>]  - g [<=°]  a 

-i  . (52) 

&2  M - g [°°]  2 

Only  the  a-g[e]  dependence  in  k,  therefore,  yields  a nonlinear  rela- 
tion between  the  differences  a-g[e]  at  large  strain  in  tensile  tests  with 
various  strailn  rates.  Such  nonlinear  relations  are  usually  observed  in  real 
tensile  tests.  Figures  7a,  7b,  9 and  10  illustrate  this  property. 

The  limits  (51)  and  (13)  are  valid  only  at  e -♦  <*>.  However,  the  stress 
ratio  and  the  stress  slope  are  rapidly  asymptotic  to  these  limiting  values 
[(51)  and  (18)],  as  shown  in  Figs. 2,  5 through  10.  Moreover,  the(a-g[e])- 
dependence  appears  to  enhance  the  rate  with  which  these  limits  are  obtained. 

From  (24)  and  (25)  we  obtain  in  the  case  of  an  instantaneous  change  of 
strain  rate 

o [°°]  - g [«°] 

— = (E-g'M)ff.  (53) 

k [cta  [®]  - g [co] ) 2 

and 

a [°°]  - g [®] 

= (E  - g 7 [“]  )ar  . (54) 

k[a  [«.]  -g[»]] 


From  our  previous  r strictions  upon  the  function  k we  have  a = a in  this 

A B 

case. 

• , do"  /v  ^ * 

Aside  from  insuring  that  — = E for  cr-g[e]  =0,  setting  g [e]  =E  has  the 

following  consequences  for  the  test  involving  a change  in  the  strain  rate, 

see  Fig. 3.  Suppose  the  function  g[e]  has  the  shape  of  a stress-strain  diagram 

/N 

da  / 

as  shown  in  Fig.  1,  then  we  know  from  (16)  and  J[18 ) that  — approaches  g [e] 

daB 

in  a constant  strain  (stress)  rate  test.  If  — <*=’  g7[e]  and  g [e]  « E then 

d€ 

If  CTj  [oo]  and  as  [®]  are  finite  then  k [8i  [°°],“]  may  not  be  equal  to 
k [oa  [<*>],<*>] . This  may  occur  if  g7!®]  = 0. 


i 


I 
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the  second  term  on  the  right-hand  side  of  (30)  is  close  to  the  value  E.  If 

the  strain  rate  is  now  increased  instantaneously  by  one  or  several  orders  of 

da 

B + 

magnitude  then  the  slope  — — (e -•  e ) in  (31)  will  be  approximately  equal  to  E 

de  o 

since  the  second  term  in  (31)  becomes  negligible.  Therefore  the  instantaneous 
change  in  slope  predicted  by  (1)  is  approximately  the  elastic  slope  E.  The 
above  sp  cialization  insures  that  this  behavior  is  independent  of  strain. 

B)  The  Derivatives  m7  = -r-p— — < 0 amd  k7  < 0 

d (a  - g [e  ] ) 

This  specialization  provides  a sufficient  condition  for 

s g7[e]  in  a tensile  loading  provided  e>0,  E=g7[e]  ;>g7[e]  and  e = a 

e[t]  i 0 in  a creep  test  provided  g 7 [e ] ^ 0 
a[t]  ^ 0 in  a relaxation  test 
de  [t] 


dcr„ 


^ 0 in  (41),  the  creep  curves  do  not  cross,  i.e., 

eft]  s e [t]  for  a.  > a„  > 0 and  every  t. 
CT1  2 1 2 


d?  ^ ^ 0 the  relaxation  curves  do  not  cross  and  a[t]  £ a[t] 

O _ . 12 

for  eL  * e2. 

Further  the  no  crossing  conditions  for  the  creep  rate  curves  2 0^ 

and  for  the  relaxation  rate  curves  s o)  appear  to  be  satisfied. 

It  can  be  demonstrated  that  for  k either  constant  or  depending  upon 

a - g [e  ] with  — 'Wr  = k7  ^ 0 and  g 7 [e]  = E ^ S g 7 [e]  provided  e = o>  0 

d (a  - g [e  ] ) ’ de 

*/  / 

as  shown  in  Figs.  2,  7a,  7b  and  10.  However,  if  the  difference  g - g is  not 

* 

always  positive  or  if  g[e]  is  nonlinear,  then  wavy  stress-strain  curves  may 

* 

result  as  shown  in  Fig. 5.  In  proceeding  from  Fig. 2 to  Fig. 5 only  g was  altered 

from  a linear  to  a nonlinear  relation.  Further  examples  of  causes  of  waviness 

are  given  in  Figs. 6 and  8.  If  k is  neither  constant  nor  dependent  on  a-g[e] 

*/ 

wavy  stress-strain  curves  may  also  occur  even  if  g is  constant  as  shown  in 
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Fig. 9.  Note  that  waviness  disappears  when  k is  made  to  depend  on  a- g [el  as 
was  done  in  Fig. 10. 

Waviness  of  the  stress-strain  curves  can  also  be  obtained  by  making  the 

equilibrium  curve  g "wavy"  as  demonstrated  in  Figs. 7a  and  7b.  In  both  cases 
* 

g = Ee.  The  g-curve  is  constructed  in  such  a way  as  to  resemble  the  yielding 
of  carbon  steel.  In  Fig. 7a  k = 3.2  s.  The  k-function  of  Fig. 7b  taken  from  [16] 
depends  on  a-gfe]  and  the  a-curves  for  various  strain  rates  follow  the  g [ e ] — 
curve  until  the  dip  disappears  at  a = 10  1 s \ Owing  to  the  constant  value 

„ -4  -1 

of  k only  the  cr-curve  for  O'  = 10  s follows  g[e]  in  Fig.  7a.  The  curves  for 

* 

the  higher  strain  rates  are  coincident  with  g = Ee.  They  will  become  parallel 
to  g[el  at  a-g[e]  values  determined  by  (19).  These  values  are  beyond  the 
limits  of  Fig. 7a. 

A comparison  of  Figs. 7a  and  7b  shows  the  importance  of  the  a-g[e]- 
dependence  in  k to  provide  for  an  early  "locking  in”  and  for  the  nonlinear 
spacing  of  the  a-curves  obtained  at  various  Oi,  see  Eq.  (51). 

Discussion 

The  model  (1)  subject  to  specializations  A)  and  B)  reproduces  in  a unified 
way  many  of  the  characteristic  features  of  metal  deformation  behavior  generally 
referred  to  as  isothermal  dynamic  plasticity  or  viscoplasticity.  These  include 

• Uniform  slope  at  the  origin 

• Linear  elastic  or  linear  viscoelastic  behavior  at  small  strain  (stress) 

• Initial  elastic  response  upon  a jump  in  strain  rate 

• Nonlinear  spacing  of  the  stress-strain  curves  at  various  strain  rates 

• Elastic  response  for  very  fast  loading 

• Existence  of  an  equilibrium  stress-strain  curve 

• 


Creep  and  relaxation. 
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In  addition  to  the  qualitative  arguments  discussed  herein  we  demonstrate  in  [16] 
by  numerical  experiments  which  simulate  constant  rate  loadings  with  various 
stress  (strain)  rates,  with  tension- tension  cycling  and  creep  and  relaxation 
conditions  that  the  general  appearance  of  the  numerically  obtained  responses 
is  very  reminiscent  of  the  actual  behavior  of  metals.  We  therefore  see  that 
a nonlinear  viscoelastic  constitutive  equation  such  as  (1)  subject  to  the 
specializations'  A)  and  B)  can  reproduce  viscoplasticity  very  well.  However, 
the  correspondence  ceases  to  hold  if  cr-g[e]  would  change  sign.  At  this  point 
modifications  will  be  necessary  such  as  discussed  in  [11] . A forthcoming 
paper  will  address  this  subject. 

Making  the  functions  depend  upon  cr-g[e]  renders  (1)  in  essence  an 
overstress  model  [17, 18]  which  has  been  considered  many  times  [19] . However, 
the  previous  use  of  the  overstress  concept  is  either  tied  to  the  usual  de- 
composition of  the  strains  into  elastic  and  plastic  parts  and  the  yield  surface 
concept  [18, 19]  or  to  very  specific  functions  [17].  Here  we  have  demonstrated 
that  none  of  these  concepts  are  necessary  to  reproduce  metal  behavior  quali- 
tatively as  long  as  there  is  no  substantial  unloading  involved.  Moreover 
creep  and  relaxation  are  included  in  (1)  in  a unified  way.  Also  we  have 
explored  and  shown  the  limits  obtained  from  (1)  for  various  extreme  conditions 
which  to  the  authors'  knowledge  has  not  been  done  before. 

We  note  that  isochronous  creep  curves  which  are  a common  way  of  repre- 
senting creep  data  in  engineering  [20,21]  can  be  obtained  from  (1).  Iso- 
chronous creep  curves  are  vs.  e cross-plots  (where  is  a parameter)  of 
regular  creep  curves  for  constant  creep  time. 
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Considering  (41)  and  the  derivation  of  it  we  see  that  it  represents  for 

t = constant  the  reciprocal  of  the  slope  of  the  isochronous  creep  curve. 

Condition  (42)  ensures  that  ^ 0;  i.e.,  the  slope  of  the  isochronous 

o 

creep  curve  will  be  positive.  Isochronous  creep  curves  of  metals  [20]  and 
plastics  [21]  have  this  property. 

In  the  interpretation  of  the  limit  obtained  in  (46)  the  condition  t-°= 
refers  physically  to  a long-time  creep  curve.  The  limit  (46)  indicates  that 
the  slope  of  this  infinitely  large  time  isochronous  creep  curve  is 

da 

-rp-  - g'[e[t  = »l]  (55) 

a® 

from  which  we  infer  that  it  must  be  equal  to  the  slope  of  the  equilibrium 

stress-strain  curve.  This  inference  is  only  valid  if  the  stress  a is  such 

o 

that  e(t  = °°)  = g [C  ];  see  Eq.  (50)  and  Appendix  III. 

These  reflections  suggest  then  that  g[e]  could  be  obtained  from  a long 
term  isochronous  creep  curve,  say  for  t = 300,000  hours  for  the  case  of  metals, 
and  could  be  obtained  from  [20] . 

In  the  above  we  have  concentrated  on  the  specializations  A)  and  B)  which 
appear  to  be  appropriate  for  metals.  However,  other  properties  of  (1) 
resulting  from  different  specializations  are  of  interest  for  the  modelling 
of  other  materials. 

If  g and  g are  selected  such  that  g/-g  / < 0 then  inverse  strain  rate 

sensitivity  results  such  as  that  shown  in  Fig. 6 which  is  obtained  from  Fig. 5 
. . * 

by  interchanging  g and  g.  It  can  be  seen  that  the  "slow"  stress-strain  curves 
are  above  the  fast  ones,  reversing  the  trend  of  Fig. 5 and  the  relation  (15). 
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It  is  also  possible  to  have  "classical"  strain  rate  sensitivity  over  one 

strain  interval  and  "reverse"  strain  rate  sensitivity  over  another  interval 
* 

by  making  g and  g intersect  once.  The  resulting  stress-strain  curves  are 
shown  in  Fig. 8. 

* 

In  Figs.  2,  5-7  the  initial  slopes  of  g and  g are  made  equal.  In  these 

cases  all  the  curves  are  initially  coincident  and  we  have  initial  linear 

* / / t 

elastic  behavior.  In  Fig. 8,  g > g initially  (k  is  constant  in  this 
example),  and  therefore  a linear  viscoelastic  behavior  results. 

Figure  8 also  illustrates  that  all  curves  with  nonvanishing  strain  rate 
have  initially  the  same  slope  as  stated  in  (17). 

If  we  begin  a creep  (a  = 0)  or  a relaxation  (e  = 0)  test  on  any  stress- 
strain  curve  then  strain  (stress)  would  always  move  towards  g.  Specifically 
in  Fig. 6,  a reversal  of  the  usual  situation  would  occur,  since  according  to  (1) 
during  relaxation  the  stress  would  increase  and  the  strain  would  decrease 
during  creep.  In  Fig. 8 a mixed  behavior  would  be  obtained  with  the  usual 
pattern  as  long  as  cj  > g and  the  reverse  behavior  for  a < g.  Finally  Fig. 9 
illustrates  that  wavy  stress-strain  curves  may  also  be  obtained  by  constructing 
special  k-functions  not  included  in  specialization  A). 
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APPENDIX  I 

Fast  and  Slow  Loading  in  Stress  Control 

To  simulate  a tensile  test  under  constant  stress  rate  3 we  set 

CT  = 3t  . 

The  solution  of  (4)  subject  to  (Al)  is  designated  by 


e [a]  = e [a,-  0]  = e [t] 


a=3t 


(Al) 


(A2 ) 


and  3 is  considered  a parameter. 

Substitution  of  (Al)  into  (4)  using  (A2)  and  (12)  yields 
CT 


e la]  = + 


f ± (-4.W  )}{exp  _ i [ — dsj  dx . 

" V^rSrvii  dx  'g[e[x]]/A  3 { e [s]  m [s.  e (s]  ] J 


g[e[a])  J g [e  [x]  ] 


For  extremely  slow  loading  3 0 and  from  (A3) 


(A3) 


e (a;  0]  = 


ae [a; 0] 
g[e  [a;0] 


Consequently 


a = g [e  [ct,-0]  ] 


(A4) 


(A5 ) 


where  e[a;0]  = lim  e[a;31.  Similarly,  we  define  e [ct  ;<*>]  = lim  e[o;3l  and 
3-0  3-“ 


proceeding  formally  we  obtain 


e fcr 


- - J ^ 


dx 


o g [s  ;°°]  ] 


and  therefore 


a = g [e  [cr[°°]  ] ] . 


(A6 ) 


(A7) 


Comparison  with  the  results  for  strain  control  shows  that  g(e]  and  gfe] 
are  the  response  curves  for  infinitely  slow  and  infinitely  fast  loading, 
respectively,  irrespective  of  the  type  of  loading. 


A2 


APPENDIX  II 


We  want  to  emphasize  that  (A9)  is  valid  only  if  the  values  of  the  functions 
at  infinity  exist  and  are  finite. 


A3. 


APPENDIX  III 


Creep  Behavior 

We  subject  (4)  to  the  following  stress  history 


cj=CT[t]  0 £ t < t with  o[t  ] = a 

o o o 


a = 0 


t £ t 


where  a is  some  smooth  function  of  time, 
given  by 

t 

>[t]  _ ct[t] 


The  total  strain  for  t ^ t is 

o 


e[t] 


-I(f7 


o g [e [t]  ] 


m[CT[T],e  [t] 


tK-(J  «°> 


which  may  be  written  as 


e [t]  = exp 


f 9t£  ts]  ]ds  J 

J e [s]m[o  e [s]  ] l o 
t ° 

o t tQ  . 

P ( P g[e[s]]ds  \ dT 

+ °o  J Vexp  “ J e [s]m[a  e [sTT/  m[a  e[T)} 

*-  rr  O SJ  J 


(All) 


Performing  the  limit  t -*  00  for  (All)  using  the  methods  of  Appendix  II  shows 

e [®]o 


e e 111  - iTTHT 


(A12) 


so  that 


e[°°]  = g 1[dQ]  . (A13 ) 

g-inverse  in  (A13)  may  not  exist.  In  this  case  lim  e = » and  creep  does  not 

t-“> 

terminate.  (Such  a case  would  be  represented  by  the  condition  g'(e)  = 0 for 
e > e and  Oq  > g[e].)  In  the  case  where  the  inverse  of  g in  (A13)  exists, 
creep  terminates  at  e = g and  we  have  only  primary  creep. 


A4 . 


It  is  easily  seen  from  (All)  that  e (t)  ^ 0 for  positive  e where  equality 
applies  only  when  ctq  = g [e [»] ] . 

The  above  concerns  a constant  stress  creep  test.  We  will  now  consider 
a constant  load  creep  test.  Then  we  enforce  a continuous  stress  history  with 
a possible  jump  discontinuity  in  O'  at  t = t 


o = o(t) 


0 £ t < t 


with  o(t  ) = ct 
o o 


for  t ^ t 


o a a A 

_ 00  , • o o • 

a = — - — and  ct  = - — — A 


where  A represents  the  cross  section  (which  is  in  general  a function  of  time 

through  e)  of  the  bar  at  time  t and  A = A[t  ] . 

o o 

Substitution  of  the  above  into  (4)  yields  for  t S t 


e [t]  = exp  - J 


g [e  [s]  Ids 


CT  A 


to  e [s]m[^~-°  , e [s] 


{ 


a*  J i(. — i 

O O J A \ rCT  A 


Vremi 


■)(.xp  - J 


If  A [t] 
is  obtained. 


and  A[t]  are  finite  and  A[o°]  # 0,  the  limit 

Proceeding  formally,  with  e [°°]  = lim  e[t], 

t“*°° 


g [e  [s]  ]ds 

pCT  A 

e [s]m|  -°A  ° , e [s] 


- ))dT 
J 

(A14) 


for  large  times  of  (A14) 
we  have 


CT  A 
o o 

A [°°] 


g[e  [°°]  1 


■G 


CT  A 
O O 


[00] 


>- 


e [CD]  ^ [CO] 


A [®]  g 7 [e  [“]  ] 


(A15) 


I 


A5. 


Again  there  may  not  be  a finite  e [<*>]  which  satisfies  (A15).  In  this  case  the 
creep  strain  becomes  infinite.  If,  however,  a finite  value  of  the  strain 
exists  then  for  a given  a at  t = 00  the  constant  load  creep  strain  is  greater 
than  or  equal  to  the  constant  stress  creep  strain  since  A < 0 and  g'  > 0. 
Physically  this  means  that  the  strain  in  a constant  load  creep  test  is  larger 
than  the  one  in  a constant  stress  creep  test. 

Equations  (34),  (All)  and  (A14)  show  a dependence  on  the  difference 
o-gfe]  attained  at  the  start  of  the  relaxation  or  creep  test.  This  difference 
is  influenced  by  the  loading  rate  up  to  the  strain  (stress)  level  at  which  the 
relaxation  (creep)  test  starts.  This  difference  is  then  multiplied  by  a factor 
with  exponential  decay  in  time.  It  is  therefore  clear  that  the  evolution  in 
time  of  the  theoretical  relaxation  (creep)  curves  is  influenced  by  the  loading 
rate  up  to  the  strain  (stress)  level  of  the  test.  Experiments  show  such  in- 
fluences [23]  . 

It  is  important  to  note  that  both  the  constant  stress  and  constant  load 
creep  test  involve  an  additional  integral  term  compared  to  the  relaxation 
case,  see  Eqs.  (34)  and  (All),  (A14). 

The  case  of  secondary  and  tertiary  creep  will  be  examined  in  detail  in 


the  future. 
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Figure  2 


Figure  3 


Figure  4 


Figure  5 


Figure  6 


Figure  7a 


Figure  7b 


Figure  8 


Schematic  Showing  the  Solution  Properties  of  Eq.  (13).  For 
given  gfe]  and  g[e]  the  curve  for  all  strain  rates  are  between 
these  two  curves. 

Integration  of  (13).  For  convenience  k = const  = 16  (units  of 
time).  Strain  rates  have  the  reciprocal  time  dimension  of  k. 
The  behavior  predicted  in  (18)  and  (19)  is  obtained  for  low 
strain  rates.  Owing  to  the  constant  k large  values  of  a-g[e] 
are  necessary  at  high  strain  rates  before  the  limit  is  reached. 
The  S'  [ e ; O']  curves  confirm  Eqs.  (14),  (15)  and  (18). 

Schematic  Showing  the  Response  to  a Change  in  Strain  Rate. 

At  A the  strain  rate  is  changed  from  0^  to  a 2 with  or ^ > or^. 

Note  that  we  distinguish  between  the  a - e and  the  CT-time 
responses. 

The  Influence  of  the  cr-Dependence  of  the  k-Function  on  the 
Relaxation  Behavior  at  eQ  = .6%.  Curve  1 is  for  linear 
relaxation  k = 16 . In  curve  2 the  k-function  used  in  Fig. 10 
is  employed  to  yield  nonlinear  relaxation. 

Stress-strain  diagrams  at  Various  Strain  Rates  Measured  in 
Reciprocal  Time  Units  of  k.  Same  conditions  as  in  Fig. 2 
except  for  a change  in  g.  Wavy  stress-strain  curves  result  for 
|ome  strain  rates  because  of  the  particular  combination  of 
g and  k. 


Same  as  Fig. 5 Except  the  Representation  of  g and  g have  been 
Interchanged,  such  that  g'  - q'  £ 0.  As  a consequence  the  curves 
for  low  strain  rate  are  above  those  for  high  strain  rate 
simulating  inverse  strain  rate  sensitivity. 

Stress-strain  Behavior  for  k = 3.2  s and  a g-curve  with  Dip. 

Note  that  the  curves  for  strain  rates  of  lO-1'  and  10°  s-1  follow 
the  straight  line  of  the  g-curve  within  the  limits  of  the  graph. 

Same  as  Fig. 6 a Except  that  k in  Units  of  Seconds  and  Taken 
from  [16]  is  Made  to  Depend  on  cj-g[e].  The  nonlinear  spacing 
of  the  5-curves  at  various  strain  rates  [in  s-1  ] are  obvious, 
see  Eq.  (51).  Note  that  the  dip  disappears  at  high  strain  rates. 

¥r 

The  Curves  g and  g Cross  and  g7-g  ^ 0 Initially  but  Changes 
Sign  at  Some  Strain,  after  which  g/-g/  ^ 0.  Therefore  a 
combination  of  normal  and  inverse  strain-rate  sensitivity  is 
obtained.  Although  g'^0)  # g ^ (0 ) all  the  curves  depart  from 
the  origin  with  slope  g7(0),  Eq.  (17). 
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Figure  9 


Figure  10 


A Particular  Limiting  Choice  of  k so  that  k[“]  = 0, 
Consequently  3-g[e]  - 0 for  all  Strain  Rates.  Although 
g = E = const,  wavy  stress-strain  curves  result. 

Same  as  Fig.  9 Except  k is  Made  to  Depend  on  a-g[e]  without 
Changing  the  Form  of  the  Function.  The  waviness  disappears 
(see  specialization  B)  and  the  curves  "lock  in". 
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